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1 Introduction 

The conformally flat perfect fluid solutions of Einstein equations with nonzero 
expansion were considered by Stephani 1 and are usually called Stephani Uni- 
verses. These solutions were rediscovered by Barnes 2 who obtained them as 
the conformally flat class of irrotational and shear-free perfect fluid spacetimes 
with nonzero expansion. 

In order to generalize the cosmological principle, Bona and Coll 3 looked 
for spacetimes admitting an iso-invariant synchronization, and without any 
hypothesis on the energy tensor, they likewise found the Stephani Universes. 

An iso-invariant synchronization implies that an irrotational and shear- 
free observer u with non zero expansion exists, and the induced metric on 
the orthogonal hypersurfaces are of constant curvature. Then, there exists an 
adapted coordinate system, u = ^9 t , such that the line element writes: 

ds 2 = -a 2 dt 2 + fPSijdxtdxi 

^ ^_R(t) a ^Rd R ]nQ C 1 ) 

1 + 2b(t) ■ f+ \K{t)r 2 

R(t), b(t) and K(t) being 5 arbitrary functions of time. 

The metric (1) is a perfect fluid solution, T = (p + p)u ® u + pg, with 
energy density and pressure given by 

iR 2 3 , T , 3 dno 

P = lP+iT^ K - Ab ^ p = - p -R~f (2) 
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Moreover, the expansion of u and the curvature of the spatial metric are ho- 
mogeneous, and they are given, respectively, by 

0(t) = ^O, K (t) = ±- 2 (K-^) (3) 

Stephani Universes are perfect fluid solutions of Einstein equations, but 
what is the physical meaning of these perfect fluids? Bona and Coll 3 gave 
the first step to answer this question: they studied when these perfect fluid 
solutions admit a thermodynamic scheme, and they showed that the ther- 
modynamic Stephani Universes admit at least a three-dimensional group of 
isometries. This topic has been again considered later 4 and, in a recent paper, 
Sussman 5 analyzes a family of spherically symmetric Stephani Universes that 
may be interpreted as either a classical monoatomic ideal gas or a matter- 
radiation mixture. 

Here we determine all Stephani Universes that represent the evolution in 
local thermal equilibrium of a (generic) ideal gas. In section 4 we present the 
main results that are based on previuos ones: firstly, on the above cited work 
by Bona and Coll 3 about the thermodynamic schemes in this class of solutions, 
that we shorten here in section 2; secondly, on our hydrodynamic character- 
ization of an ideal gas in local thermal equilibrium 6 , that we summarize in 
section 3. 

2 Thermodynamic Stephani Universes 

The usual physical interpretation of a perfect fluid follows when it has a con- 
servative evolution in local thermal equilibrium (l.t.e.). This means: 

• Energy-momentum conservation: V • T = 0. 

• The energy density p is decomposed in terms of the matter density r and 
the specific internal energy e: p = r(l + e). 

• The equation of conservation of matter is required: V • (ru) = 0. 

• The thermodynamic variables are related by equations of state compat- 
ibles with the thermodynamic relation: Tds = de + pd(l/r). 

Bona and Coll 3 looked for the Stephani Universes admitting the above ther- 
modinamic scheme. The results that will be uscfull for us here can be summa- 
rized in the following two lemmas. 
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Lemma 1 A Stephani Universe admits a barotropic thermodynamic scheme 
iff it admits a six- dimensional isometry group. 

So, this case corresponds to the Friedmann- Robertson- Walker limit that oc- 
curs when K = constant and b = constant. 



Lemma 2 A Stephani Universe admits a strict (non barotropic) thermody- 
namic scheme iff it admits a three-dimensional isometry group. 

For the thermodynamic Stephani Universes, the metric may be written 

ds 2 = -c?dt 2 + n 2 8i dx i dx j 

n = -L, a = Rd R InL (4) 

z 

= R(t) = z 

l + K^w 1 ^ ~ 1 + |r 2 

R(t) and K(t) being two arbitrary functions of time. 

The symmetry group is spherical, plane or pseudospherical depending on e to 
be 1, or —1. 

For the thermodynamic perfect fluid solutions (4) the expressions (2) for 
the energy density and pressure become 

3R 2 3 . T o. 3 Oro ,_. 

P= lP + R-^- K ^ V=- p -R~i (5) 



3 Ideal gas evolving in local thermal equilibrium 

The definition of local thermal equilibrium given at the begining of the previ- 
ous section introduces thermodynamic variables, like r, e, s, and T, that are 
not present in the perfect fluid energy tensor. Then, a natural question arises: 
does an equivalent formulation of the l.t.e. exist that reduces to the addition 
to the energy-momentum conservation equation of a new condition involving 
only the hydrodynamic variables (u,p,p)7. We gave years ago 7 a positive an- 
swer to this question: 

Lemma 3 A divergence-free perfect fluid energy tensor evolves in l.t.e. if, and 
only if, the space-time function x = p/p, called the indicatrix of l.t.e., depends 
only on the variables p and p : d\ A dp A dp = 

This result has the conceptual interest of giving an exclusively hydrody- 
namic version of l.t.e.. But moreover it also has a practical utility because 
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provides a deductive algorithm to detect wether or not a given divergence-free 
perfect fluid energy tensor evolves in l.t.c.. In the Spanish Relativity Meeting- 
96 8 we pointed out these applications and gave a concrete exemple studying 
the thermodynamic class II Szekeres-Szafron space-times. 

Nevertheless, in order to have a more accurate physical meaning of a per- 
fect fluid it will be usefull to study the indicatrix function of several particular 
media. Elsewhere 6 we have given this hydrodynamic characterization for some 
representative cases. We present below in two propositions a part of our results 
about ideal gases that we will use in next section. 

A (generic) ideal gas satisfies the equation of state p — krT. For it we 
have the following hydrodynamic characterization 6 : 

Proposition 1 The necessary and sufficient condition for a non barotropic 
and non isoenergetic (p ^ 0) divergence-free perfect fluid energy tensor (u, p,p) 
to represent the l.t.e. evolution of an ideal gas is that the indicatrix function 
X = p/p be a function of the variable n = p/p, x = xM 7r; 

d-X A dir = 0, x ^ 7T (6) 

Proposition 2 A non barotropic and non isoenergetic divergence-free perfect 
fluid energy tensor (u,p,p) verifying (6) represents the l.t.e. evolution of the 
ideal gas with specific internal energy e, temperature T, matter density r, and 
specific entropy s given by 

e(p,p) = e(7r)=e(7r)-l; T(p,p) = T(tt) = £e(7r) (7) 

r (p,P) = -rr; s{p,p) = s + k[l (j){-K)d-K -\np] (8) 

e(7r) J 

e(-7r) and (j)(7r) being, respectively, 

e(w)=e efo*W d ", ^ ^ — — ^— — -; = — i (9) 

(x(tt) -7r)(7r + l) XW-tt 

The above two propositions provide an algorithm to detect, in any given family 
of divergence-free perfect fluids T = {T = [u a (x^), p(x l3 ),p(x^ l )]} 7 those that 
represent an ideal gas evolving in local thermal equilibrium: 

1. Calculate the coordinate dependence of the space-time functions p/p = 

x(x^) and p/p~Tr(x^) 

2. Determine the gas ideal subset of T by imposing the hydrodynamic con- 
dition of the proposition 1: d\ A dn = 0. 
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3. Obtain the explicit expression of the indicatrix function: \ = x( 7r )- 

4. Calculate, from x(7r), the functions e(7r) and VK 71 ") given in (9), and obtain 
the thermodynamic variables using (7) and (8). 



4 Ideal gas Stephani solutions 

In order to determine the ideal gas Stephani Universes we start from the ther- 
modynamic ones (4), and we must study the restrictions that our ideal gas 
characterization imposes on the functions R(t) and K(t). We follow the algo- 
rithm presented at the end of previous section: 



1. Variables \ an d ^ 

From (5) and (4), a direct calculation lends to: 

p/p - x(tt, J?) = 7T + - - -(tt + + - + (tt + 1 - a)-^—] (11) 

o o (X (X Qj J\ 

where the prime indicates derivative with respect to the variable R. 

2. Ideal gas hydrodynamic condition: d\ A dn = 

This condition restricts, in a first step, the functions p(R) and K(R) that turn 
out to be submitted to the second order differential equations 

// 2 «' 2 , K " 
a"~ — =a'- (12) 

K"'-^-\) = ™"- a (13) 

Secondly, we put every solution p(R), K{R) for these equations in the ex- 
pression (5) of the energy density p and we obtain a generalized Friedmann 
equation for the expansion factor R(t): 

p(R) = ^- + ^- K2 m (14) 
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3. Indicatrix function: \ = x(7r) 

At this point we look for the functional expression to the indicatrix function. 
Assuming equations (12) and (13), the indicatrix (11) becomes 

x(tt) = an 2 + 7?r + (3, a + ft = 7 - 2/3 (15) 

4- Thermodynamic variables 

Finally, the form (15) of the indicatrix function determines, using the results 
in proposition 2, the other thermodynamic variables, and so, the properties of 
the associated thermodynamic scheme. 

In contrast to the Sussman 5 intuitive approach leading to very partial 
results, our algorithm gives directly, and in the above simple form, the whole 
class of ideal gas Stephani Universes. After that, we may look for a more 
specific ideal gases by imposing the corresponding restrictions to the indicatrix 
function (15). For example, if one collapses the whole space of solutions of 
equations (12) (13) to the particular solutions a(R) = Constant, and one 
imposes the gas to be monoatomic, one obtains the partial results by Sussmann. 
The physical meaning of the other cases will be considered elsewhere. 
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